1 



Is there elliptic flow without transverse flow? 
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^ ■ Azimuthal anisotropy of final particle distributions was originally introduced as a signa- 
O , ture of transverse collective flow We show that finite anisotropy in momentum space 
Q . can result solely from the shape of the particle emitting source. However, by comparing 
'nI" I the differential anisotropy to recent data from STAR collaboration we can exclude 
^ ' such a scenario, but instead show that the data favour strong flow as resulting from a 
hydrodynamical evolution [0,0. 

1. Bjorken Brick 

Q . 

2 ■ The intensity of particles emitted from a steady thermal source depends on the size 
of the source facing the direction of emitted particles. If the shape of such a source 
is not azimuthally symmetric, the particle yield is anisotropic too. Thus it is easy to 
^ ' parametrize a source with an essentially arbitrary value of the elliptic flow coefficient V2, 
' but the transverse momentum dependence of the differential V2{pt) is much less trivial. 

To gain insight what kind of elliptic flow might result from surface effects only, we 
constructed a simple parametrization of the source. Let us assume that the longitudinal 
expansion of the source is boost invariant, but that the source does not expand transver- 
sally at all. To simplify the calculation we take the source to be rectangular in transverse 
plane and to have a constant temperature T. Particle emission now takes place on the 
surfaces of the system during its lifetime (from thermalization time tq until time t) and 
from its entire volume at flnal breakup at time t (see Figs. |^ and Unlike in where 
the source was parametrized using a Gaussian proflle, our source contains emission also 
through surfaces which normal vectors are spacelike. If these surfaces are opaque, this 
allows us to obtain flnite V2 even if transverse flow is absent. 

We calculate particle distributions using a modifled Cooper- Frye procedure 0: 

= / • T) e{da,pn, (1) 

where the 0-function is needed to remove negative contributions to particle number. For 
the freeze-out surface described above, the integral can be calculated analytically. Using 
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Figure 1. Boost-invariant expansion of the Figure 2. Projection of the Bjorken brick 
source and emission in two stages. to xyt-speice. 



the Boltzmann approximation and the customary definition of the coefficient V2, we get 

V2(Pt) = 7 7 \ ^ 7 ^ T- 2 

3 {7itmtK,{^) IJy + (t2 - r2)piii'o(ff) (/, + 

Differential elhptic flow is now essentially controlled by three parameters: the ratio of 
the system size in x- and y-directions Ix/ly, the ratio of lifetime to spatial dimension t/ly 
and temperature T. The effect of thermalization time Tq can be compensated by changing 
the lifetime and therefore we prefer to keep it fixed to value tq = 0.6 fm/c inspired by our 
hydrodynamical calculations 0]. 

In Fig. ^a) the differential elliptic flow of pions obtained from this model using four 
different parameter sets is shown. The solid line corresponds to a case where parame- 
ter values are motivated by the values from hydrodynamical simulations of semi-central 
collisions {Ix/ly = 0.58, t/ly = 1.35, T = 140 MeV). The qualitative behaviour of V2 
is surprisingly similar to that seen in the Low Density Limit approximation where 
the anisotropy is estimated by calculating a flrst order correction to the free streaming 
distribution. At low values of transverse momenta, f 2 increases rapidly but saturates al- 
ready around pt = 200 - 300 MeV/c. As demonstrated by the dashed line obtained using 
Ix/ly = 0.5 and t/ly = 1, a change in the ratios scales the saturation value of V2, but does 
not affect the value of pt where saturation takes place. 

The only variable which affects the slope of V2{pt) is the temperature of the system 
as shown by the dotted line (T = 500 MeV). The effect is weak and even outrageous 
values of parameters {Ix/ly = 1/15, T=75 GeV, dashed-dotted line) fail to reproduce the 
data. Therefore we may conclude that the present data cannot be described without the 
presence of strong transverse flow at RHIC. 

The generalization of the Bjorken brick to allow transverse flow is straightforward al- 
though the resulting expression for V2{pt) is much more complicated than Eq. (0). In addi- 
tion to the parameters mentioned earlier, we have the flow velocities in x- and y-direction 
Vx and Vy. This generalized parametrization allows us to check which one contributes 
more to the observed anisotropy: anisotropic flow fleld or source shape. Figure ^b) shows 
our fits to the data assuming that both source shape and velocity field are anisotropic 
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Figure 3. Differential elliptic ffow of pions from Bjorken brick without flow (left) or with 
additional flow (right) compared to elliptic flow of charged particles at RHIC 0. For 
explanation of the curves, see the text. 



(solid line, = 0.6, Vy = 0.57, Ix/iy = 0.93), only the velocity field is anisotropic (dashed 
line, Vx = 0.6, Vy = 0.57, Ix/ly = 1) or only source shape is anisotropic (dotted line, 
Vx = Vy = 0.7, Ix/ly = 0.73). Due to the crudeness of the model the differences between 
these fits are not meaningful, but the data allows all three different possibilities. However, 
an isotropic velocity field from a deformed source depicts very interesting behaviour of 
saturating V2 around pt = 1.5 - 2 GeV/c, in the same region where data do. Whether this 
is anything but a coincidence and whether it is possible to produce this kind of behaviour 
using a dynamical model remains to be seen. 



2. Hydrodynamic calculations 

We have used a boost invariant hydrodynamical model ^ to describe the time evolution 
of the collision system at RHIC P,|5|JT0[]. In our approach the initial conditions of the 
hydrodynamical evolution were fixed by scaling the initial densities used to reproduce 
SPS data ^ until the final particle multiplicity measured by PHOBOS collaboration ^ 
was reached. In the most central collisions this leads to maximum energy densities of 
e = 22.3 - 23 GeV/fm^ and temperatures of T = 270 - 330 MeV, depending on the 
equation of state (EoS). 

In Figs. ^ and § we show our results for minimum bias collisions obtained using two dif- 
ferent EoSs — one with a phase transition to QGP (EoS Q) and one without (EoS H). The 
agreement between the charged particle data and hydrodynamic calculation is excellent 
up to Pt ~ 1.5 GeV/c where the measured elliptic flow begins to show signs of saturation 
to some value whereas the hydrodynamic result keeps increasing. The pion data show 
similar agreement independent of the EoS while at first glance the proton V2{pt) data seem 
to favour the existence of a phase transition. However, as shown in [p|,p!0|, proton elliptic 



flow is sensitive to freeze-out temperature and details of initial distributions, contrary to 
the elliptic flow of pions and charged particles. Therefore no conclusions about the EoS 
can be drawn without further constraints on the model. The data also show a strong 
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mass dependence of V2{pt) which is typical to hydrodynamical models: the heavier the 
particle, the smaller V2 at small values of pt- 

As a function of centrality (see our fit to data is less successful. As can be expected, 
we can reproduce the data in central and semi-central collisions, but we overestimate 
elliptic flow in peripheral collisions, due to the loss of a sufficient degree of thermalization. 
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Figure 4. Differential elliptic flow of Figure 5. Differential elliptic flow of pions 
charged particles in minimum bias colli- and protons-|-antiprotons in minimum bias 



sions at RHIC compared to data p[. 



collisions at RHIC compared to data 0. 



The observed elliptic flow depends on the amount of rescatterings during the expansion. 
Hydrodynamics assumes zero mean free path, or equivalently, infinitely strong rescatter- 
ing. It therefore provides practical upper limits for the anisotropy of particle distribution. 
That the experimental data reach these limiting values for elliptic fiow and also refiect the 
mass dependence typical for a hydrodynamic system points towards early thermalization 
and a subsequent hydrodynamic evolution in central and semi-central collisions at RHIC. 
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